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RESUMEN 

Se discute en este articulo la forma que adoptan los espectros para diseno en la actual 
reglamentacion para diseno sismico en la ciudad de Medellin. En particular, se discute y proponen 
nuevas formas para la zona descendente de los espectros para diseno y control de danos, que 
adolecen de las mismos problemas que los espectros que se proponen en la actual norma 
colombiana NSR-98 y en las reglamentaciones actuales de muchos otros paises del mundo. 

En la zona descendente de los espectros para diseno, la reglamentacion actual sobrestima 
sistematicamente y de manera creciente con el periodo estructural el desplazamiento maximo 
esperado. Mas grave todavia cs que las normativas en discusion violan los principios de la 
dinamica estructural cuando postulan desplazamientos maximos del terreno infinitos. 



ABSTRACT 

The form that adopt spectra for seismic design, under current regulations in the city of Medellin, 
are discussed in this article. In particular, new forms are proposed and discussed for the descendent 
zone of service and design spectra, which suffer the same problems that spectra proposed under 
current NSR-98 Colombian regulation, and others from many other countries in the world. 

In the current regulation, the descendent zone of the seismic spectra systematically exceed and, 
in an increasing way with the structural period, the maximum displacement expected. Even more 
serious is that current norms violate the principles of structural dynamics when they postulate 
infinite maxima displacements of land. 
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INTRODUCCION 

La norma colombiana NSR-98 para diseno si'smico de edificaciones, como muchas otras 
reglamentaciones alrededor del mundo, especifican el sismo para diseno a traves de un espectro de 
diseno de seudoaceleraciones. En el caso mas general, este espectro. por simplicidad en su 
formulacion, se ha dividido en tres partes: la zona de periodos bajos, definida a tra%e$ de una recta 
que parte de la aceleracion maxima esperada en la superficie del suelo y termina en la aceleracion 
espectral maxima; la zona de periodos intermedios o periodos de amplificacion maxima. der”r.:da 
por una constante que representa la aceleracion espectral maxima especificada . La tercera rama. :■ 
rama descendente, tradicionalmente se ha especificado proporcional a IT'. En el Codiz: 
Coloinbiano de Construcciones Sismo-resistentes de 1984 se propuso para a un \alor de 2 3. que se 
cambio en la version de 1998 de esta reglamentacion a un valor de 1.0. Se podria decir que esta 
ultima version es un poco menos conservadora que la anterior, pues se propone una eaida rr.ua 
pronunciada del espectro de aceleraciones. 

Por definicion, los espectros de aceleraciones para diseno no representan el espectro de respcesia 
de aceleraciones de una serie de tiempo que pueda considerarse como el acelerograma de irser : 
mas bien se parece mas, o se fundamenta en, lo que se ha llamado espectro de amenaza ur-rrcr-rt. 
que se define como el trazo que une las aceleraciones espectrales asociadas independientemrrre a 
cada periodo estructural con una probabilidad de excedencia definida, en un periodo de uemp; 
tambien definido. En otras palabras, el trazo que une las aceleraciones espectrales asccucua u 
mismo periodo de retomo, trabajando cada periodo estructural independientemente. Es:e e-szeezr: 
de amenaza uniforme se modifica un poco para simplificar sus trazos y para considerar err i_r_r :>:s 
casos incertidumbres provenientes del calculo del periodo estructural de la edificacion. 

De la definicion de espectros de diseno se desprende que, como un todo. no mm-roc cacu 
periodo estructural independientemente, la energia total involucrada bajo ese trazo scjpen :x 
creces la energia total, como medida de la intensidad del sismo, que se calcularia supxasnac u 
misma probabilidad de excedencia y en el mismo periodo de tiempo usados para el c-tcrar ze 
espectro de amenaza uniforme. En otras palabras, el espectro de amenaza uniforme se zxnoae *er 
como la envolvente de los espectros producidos por los sismos “maximos creibles" de zaca mu zt 
las sismofuentes que afectan el sitio en estudio. Esta ultima definicion debe tomarse rm. jwi .ice 
como ilustrativa, pues rigurosamente no es claro ni el procedimiento para efectuar d z»rm ): 
descrito ni la definicion de “sismo maximo cretble”. 



Si cada edificacion se pudiera asociar estrictamente a un sistema de un solo grade de 
como la combinacion de varios de ellos, la sobre-estimacion descrita en los espeztrx 
diseno estamos utilizando no tendria ninguna incidencia, y es sobre este supuesto qor j 
metodologfa descrita para la obtencion de los espectros de diseno. En general, auoame 
sistema se representa como la combinacion de varios sistemas de un solo grado de libemd. 
ellos, llamado el fundamental, participa en una proporcion mucho mayor que los 
respuesta final del sistema. Es claro que mientras mas alta sea la participate* 
fundamental, mas justifiable es el procedimiento descrito para la obtencion del espectro 
En los casos en que esta premisa no se cumpla, se obtiene un sobre diseno que, hasta x« 
es sano, dado que estos casos normalmente ocurren en edificaciones con asimetri** a 



n: 



It 

n 




1 En algunas reglamentaciones, como es el caso de Colombia en su espectro de diseno. por 
zona plana se inicia desde el periodo nulo (T=0.0). 
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elevacion, que generan concentraciones de esfuerzos y otras medidas de la respuesta estructural en 
algunas zonas, que no estan completamente previstas en las metodologias que para diseno se usan. 

Despues de esta breve descripcion de los espectros de diseno y la justification del procedimiento 
que actualmente se sigue para estimarlos, volvamos un poco sobre las formas simplificadas 
propuestas y lo que estas significan a la luz de la Dinamica Estructural. De acuerdo con las formas 
descritas, el espectro de diseno queda completamente definido en terminos de cinco parametros: la 
aceleracion maxima a nivel de la superficie del suelo, A niax , la aceleracion espectral maxima, Sa inax , 
el primer periodo de quiebre del espectro, T 1; que define la zona de periodos bajos, el segundo 
periodo de quiebre del espectro, T 2 , que define el final de la meseta de amplification maxima, y el 
exponente, a, que define la forma de la caida del espectro. En la figura 1 se muestra esta forma 
espectral y los parametros que la definen. 




Figura 1 Forma del espectro de aceleraciones para diseno 



El argumento que ha impuesto la definition del parametro a, que controla la caida del espectro, 
parece estar asociado a la mayor seguridad que se le debe dar a las estructuras de periodos altos, que 
generalmente representan las grandes obras de la ingenierfa como son los edificios altos, puentes, 
etc. 

En cuanto a las caidas que presentarfan los espectros de amenaza uniforme, que serfa lo 
estrictamente correcto hacer, se considera que hasta el momento en el pais no hay leyes de 
atenuacion confiables en la zona de periodos altos. Para conseguirlas es necesario filtrar y corregir 
cuidadosamente los acelerogramas, como se expone mas adelante, antes de hacer regresiones sobre 
las ordenadas de los espectros de respuesta de estos registros. La otra altemativa es a traves del uso 
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de leyes de atenuacion basadas en modelos fisicos de atenuacion de las ondas y en funcion de unos 
cuantos parametros. Este ultimo tipo de leyes de atenuacion esta menos expuesto a la intromision 
de componentes de frecuencia espurias producto de la electronica de los registradores, lo que las 
hace mas apropiadas para la estimacion de espectros de amenaza uniforme en la zona de periodos 
altos. 

Aunque ya se menciono que los espectros para diseno no son exactamente iguales a los espectros 
de respuesta de acelerogramas de diseno, las formas espectrales calculadas de acelerogramas 
registrados, cuidadosamente filtrados, se asemejan a las formas espectrales de diseno, con trazcs 
mas suaves y envolventes, estos ultimos, que permiten decidir sobre los rangos de valores de 
algunos de los parametros que definen los espectros de diseno. Por lo tanto, es sensato calcular a. 
ajustando los espectros de diseno a formas espectrales propuestas de espectros de respuesta 
calculados de sismos registrados 1 . A1 hacerlo resulta claro que los valores hasta ahora propuestos 
nada tienen que ver con lo que se presenta en la realidad. Las caidas de los espectros de respuesta 
de sismos registrados son en todos los casos mucho mas pronunciadas que lo que hasta ahora se ha 
propuesto en el reglamento nacional, por ejemplo. No se ve, entonces, ningun otro argument: 
diferente al de la seguridad de estructuras de periodos altos, mencionado antes, que justitique .:s 
valores hasta ahora propuestos para el parametro a que define la caida de los espectros de diseno 

El interes repentino por tratar de encontrar un mejor ajuste al parametro a que define la caida de. 
espectro de diseno no es caprichoso: obedece a razones de mas fondo que la simple meticulesjcrbr 
en la forma espectral propuesta para diseno. En la reglamentacion de 1984 se introdujo el cccrro. 
de derivas como componente fundamental en el diseno sismo-resistente de edificaciones. oblipco: 
al calculo de los desplazamientos de la estructura. Este calculo no es novedoso en si. es r. 
resultado, que juega un papel importante en el diseno de la estructura. La ultima version re ^ 
reglamentacion para diseno sismico de edificaciones, la Norma NSR-98, tratando de contrArc =sls 
los danos a elementos no estructurales adheridos a la estructura, volvio mas estricto el centre*: re ^ 
derivas, reduciendo sustancialmente el limite permitido para estas. Los calculistas coinciden en rj-e 
actualmente, para la mayoria de las estructuras, el control de derivas es el estado limite que desire a. 
geometria y disposition de los elementos estructurales dentro de la edification. De acuerdc :i ■: 
anterior, resulta que para el diseno como lo concibe hoy la Norma NSR-98, es fundaments, ecu. 
estimacion muy precisa de los desplazamientos que tendria la estructura cuando esu se 
sometida a los sismos de diseno y control de danos. 

En la figura 2 se muestran los desplazamientos maximos que tendria un sistema de ur. grace re 
libertad cuando la respuesta estructural esta especificada en terminos del espectro de respresci re 
seudoaceleraciones (espectros de diseno) propuestos en las dos versiones que ha temdc _£ remra 
para diseno sismo-resistente en Colombia. En otras palabras, se muestran los esperrros re 
desplazamiento con base en los espectros de diseno. Para efectos de comparacion. se supree rue e, 
unico cambio que sufrieron los espectros de diseno entre las dos versiones es en el expeosne r. 
Igualmente, se incluye con un trazo punteado un espectro de desplazamiento que strre jrs 
principios de la dinamica estructural: el desplazamiento espectral es asintotico al despl^rririBSc: 
maximo del terreno, a medida que aumenta el periodo estructural. 



1 Mas adelante se presentan ejemplos de estos ajustes realizados para la Microzonitic-^.cc re 

Medellin 
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Figura 2 Espectros de diseno propuestos en las normas colombianas 

para diseno si'smico y los espectros de desplazamiento asociados 

Es claro de la figura 2, que a partir del periodo de quiebre T 2 , la reglamentacion anterior y la 
actual sobrestiman sistematicamente y de manera creciente con el periodo estructural, el 
desplazamiento maximo esperado. Mas grave todavla es que las normativas en discusion violan los 
principios de la dinamica estructural cuando postulan desplazamientos maximos del terreno 
infinitos. 

Esta discrepancia notable es evidente en toda su dimension ahora que el calculo de los 
desplazamientos es crucial en el diseno sismo-resistente de edificaciones debido a los estados h'mite 
que establece la norma. Si el diseno estuviese enteramente controlado por la resistencia, dictada por 
el espectro de seudoaceleraciones, el impacto serfa mucho menor debido a que las aceleraciones de 
diseno decrecen con el periodo, aunque se estarian igualmente sobre-estimando sistematicamente 
estas en la zona descendente del espectro. 

En lo que sigue y fundamentado en trabajos similares en otras partes del mundo y en los 
acelerogramas registrados hasta ahora por la Red Acelerografica de Medellin, se plantean nuevas 
formas espectrales, especfficamente para la zona descendente del espectro, en cada una de las zonas 
homogeneas en cuanto a respuesta slsmica que define la microzonificacion sfsmica de Medellin. 

ANTECEDENTES 

Bommer y Elnashai, 1999, por las mismas razones que impulsan este trabajo: aceleraciones 
espectrales para diseno que producen desplazamientos crecientes en funcion del periodo estructural, 
proponen espectros de desplazamiento para diseno si'smico basados en un conjunto de 
acelerogramas registrados en Europa. Su trabajo esta dirigido a la definition de leyes de atenuacion 
sobre las ordenadas del espectro de desplazamiento, que les permite deducir espectros de 
desplazamiento para diseno sfsmico. Las formas espectrales deducidas siguen los principios de la 
dinamica estructural: el desplazamiento espectral es asintotico al desplazamiento maximo del 
terreno a medida que aumenta el periodo estructural. 



6 



Juan Diego Jararmllo 



Es importante destacar del trabajo de Bommer y Elnashai el procedimiento seguido para tiltrar 
los acelerogramas de tal forma que resulten espectros confiables en la zona de periodos altos: Cada 
registro se hizo pasar por un fdtro pasa-bajos iniciando en un periodo de corte de Ins \ 
disminuyendo este, dependiendo de los resultados conseguidos, a 10s, 7s, o 5s. Para cada periodo 
de corte se realizo una doble integracion para encontrar los registros de velocidad y desplazamiento 
con el objeto de decidir, basados en la apariencia de estas trazas, la credibilidad del registro de 
aceleracion del que proveman. Cuando los registros de velocidad y desplazamiento no mejoraban 
sustancialmente en comparacion con el periodo de corte anterior, se suspendia el procesc > se 
tomaba el registro de aceleracion pasado por el ultimo filtro como el crefble para seguir el proceso 
De 183 registros de aceleracion triaxiales seleccionados para el estudio solo 123 fue posible integral 
hasta un periodo de 3.0 s. 

Tolis y Faccioli, 1999, realizaron un trabajo muy similar al de Bommer y Elnashai cce e. 
objetivo, nuevamente, de proponer espectros de desplazamiento para disefio si'smico. La piopoesta 
final resultante es muy semejante a la de Bommer y Elnashai solo que con trazos mis simr.es ;. 
suaves. 

Ordaz, Miranda y Aviles, 2000, mantienen la idea de espectros de disefio de aceleracicc. per: 
proponen una modiftcacion a la rama descendente de estos espectros, de tal manere c-c .ns 
espectros de desplazamiento o los desplazamientos maximos estructurales que se calculen a par_r 
de estos sean crefbles y se mantengan dentro de los principios de la Dinamica Estructural La :uee 
propuesta por Ordaz y colaboradores en el articulo de la referencia, con una pequena nxxt:::cac. :c. 
es la que se adopta en este trabajo. 

Igual que los autores del articulo mencionado en el parrafo anterior, el autor considera cue .. 
menos por el momento, se debe mantener la idea de espectros de disefio de acejerac: :c. 
simplemente modificandolos para ajustarlos a los requerimientos actuales de las nuevas escracegnau 
de disefio. Solo hasta que el metodo de disefio por desplazamientos este ampliamente acegcao: p:r 
la comunidad y sea adoptado como metodologfa, al menos altemativa, en las normas g-ara acse-no 
sfsmico de edificaciones, se debe proponer este tipo de espectros. El autor tambien constoera cue e. 
cambio de estrategia de disefio a lo sumo reduce las incertidumbres inherentes al process ce ucse5: 
y puede tambien ayudar a la comprension del fenomeno fisico que se enfrenta . pero s: se -esceiar. 
los principios de la ffsica, como hasta ahora se ha venido haciendo, los resultados finales no c«e-«rr 
cambiar sustancialmente. 

PROPUESTA 

Se propone en este trabajo modificar la forma de la rama descendente de los esnecmrs a; 
aceleraciones para disefio, de manera que los desplazamientos espectrales que se calculen a gamer c.e 
estos sean crefbles y se mantengan dentro de los principios de la Dinamica Estructurai 

Como se menciona en la introduccion, la rama descendente tradicionalmente se ha esnecnrscic : 
proporcional a 1 77". En el Codigo Colombiano de Construcciones Sismo-resistentes ce l-Su se 
propuso para a un valor de 2/3, que se cambio en la version de 1998 de esta reglamenuar-nn i un 
valor de 1.0. Esta forma resulta en un espectro de desplazamientos en la zona descencente ce 
espectro, proporcional a T 2 ' a , que significa que para cualquier valor de a menor que 2 se ncne un 
espectro de desplazamientos monotonamente creciente con el periodo, que viola el pnncurtc ce u_e 
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para periodos altos el desplazamiento espectral debe tender asintoticamente al desplazamiento 
maximo del terreno. 



Se propone, entonces, para la rama descendente del espectro de diseno la siguiente forma 
funcional: 



Sa = Sa„ 



k + (l-k)\ 



f t Y 







( 1 ) 



donde: 

Sa max : Aceleracion espectral maxima (ver figura 1) 

T< Segundo periodo de quiebre del espectro (ver figura 1) 

k: Parametro 

y._ Parametro 



El espectro de desplazamiento, Sd, que se deduce de esta expresion seria: 



Sd = 




SaT 2 
An 2 



(2) 



Sd = Sd„ 



k + (l-k)\ 



( T Y 

1 2 

T 

v 1 / 



( 3 ) 



donde el desplazamiento espectral maximo, Sd mcLX , es igual a: 



OJ _ S a mnx ' ^2 

max An 1 



( 4 ) 



El parametro k tiene un significado ffsico claro, es la relacion entre el desplazamiento maximo 
del terreno, D max , y el desplazamiento espectral maximo, Sd max . 



D max = Hm(Sd) 



( 5 ) 



yf. _ ^ max 



Sd„ 



( 6 ) 



El parametro ’/define la forma en que cae el espectro de desplazamiento entre el desplazamiento 
espectral maximo, Sd max , y el desplazamiento maximo del terreno, D max . Si /es igual a cero, que es 
el caso extremo, el espectro de desplazamiento es piano a partir de T 2 e igual a Sd mat . En este caso, 
el espectro de aceleraciones Sa es inversamente proporcional a T 2 , que ya es bastante diferente a lo 
que actualmente define la norma. En general, y especialmente en depositos de suelos blandos. los 
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espectros de desplazamiento muestran un pico bien definido en el periodo fundamental del deposito 
de suelo y a partir de este caen, mas fuertemente al principio y suavemente al final, al 
desplazamiento maximo del terreno, D max . Esta caida fuerte al principio esta controlada por el 
parametro y, y a medida que el periodo T aumenta se convierte en kSd max . 

La expresion que se propone para la rama descendente de los espectros de diseno (ecuacion 1 >. 
en terminos de dos parametros, y y k, en contraste con el tradicional a, permite representar las 
diferentes formas de los espectros de desplazamiento encontradas sin complicaciones adicionaks > 
sin un aumento grande de parametros. 

METODOLOGIA 

Los parametros Sa max y T 2 que definen los espectros de diseno de cada una de las micrczonas 
sismicas de la ciudad de Medellin se mantienen, es decir, no son objeto de los ajustes que se 
implementan en este trabajo. Quedan entonces por ajustar solo los parametros yy k. que denser. .a 
forma de caida en la zona descendente de los espectros de diseno. 

Como se menciono en la introduction, aunque los espectros de diseno no comestvcoer. a 
espectros de respuesta de “acelerogramas de diseno”, para tener una idea de los valores que *i:o ciz 
y y k en espectros de respuesta de registros reales, que sirvan de base para ajustar los espec~:-s de 
diseno, se estimaron estos dos parametros para todos los acelerogramas registrados hastu eli-rre rcc 
la Red Acelerografica de Medellin, RAM. 

El procedimiento adoptado en este trabajo para procesar los acelerogramas registradc-s r>:r 
RAM es enteramente similar al descrito por Bommer y Elnashai en el trabajo referenciadc Ik 
acelerogramas horizontales, solo 335 fue posible integrar satisfactoriamente luego de apikar :s 
pasa bajos con periodos de corte mayores a 5s. Este rechazo tan grande esta asociade ccc -as 
intensidades hasta ahora registradas por la red; las aceleraciones maximas registradas oscilar ercre 2 
cm/seg 2 y 10 cm/seg 2 . En estos casos, la proportion de ruido en la serial registrada es cocszcer-ane 

El parametro k se estimo en cada caso de acuerdo con la ecuacion 6, en donde el desp- i. ~ a. ~ie c: : 
maximo del terreno, D max , se leyo del registro de desplazamientos calculado. despues ce e:ecc-ar 
una doble integracion del registro de aceleraciones. El desplazamiento espectral mix ; me se e - ; 
como la ordenada maxima del espectro de desplazamientos, y se acompano esta lectura de, rer-:o: 
T 2 , en el cual se presentaba. 

El parametro yse estimo a traves de una regresion lineal de la ecuacion 3. usande Los .iL«cres ue 
k, Sd max y T 2 hallados en el paso anterior. En cada regresion se tomaron pares de vakres id. 7 • 
partir de T 2 y hasta un periodo 0.1s por debajo del periodo de corte final del registro. 
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RESULTADOS 



ESTACION 


N 


k 


ak 


Y 


cry 


EAU 


(1) 


2 


0.37 


0.05 


1.04 


0.04 


SPE 


(1) 


13 


0.33 


0.08 


1.78 


0.32 


FMI 


(1) 


12 


0.37 


0.06 


1.54 


0.63 


ECC 


(2) 


11 


0.30 


0.07 


1.43 


0.49 


MAN 


(3) 


29 


0.22 


0.05 


1.59 


0.49 


POL 


(3) 


3 


0.32 


0.09 


1.80 


0.35 


UDM 


(3) 


16 


0.34 


0.07 


1.27 


0.31 


ISJ 


(5) 


1 


0.36 


- 


1.18 




UNA 


(5) 


1 


0.26 


- 


1.81 


- 


UEA 


(6) 


14 


0.31 


0.06 


1.27 


0.38 


COM 


(8) 


3 


0.29 


0.07 


1.74 


0.12 


SOL 


(8) 


10 


0.32 


0.06 


1.20 


0.31 


EVH 


(8) 


12 


0.31 


0.06 


1.69 


0.47 


CSJ 


(10) 


12 


0.33 


0.05 


1.53 


0.35 


EMO 


(10) 


9 


0.32 


0.07 


1.66 


0.62 


EVT 


(10) 


9 


0.32 


0.06 


1.44 


0.50 


LIC 


(11) 


2 


0.31 


0.04 


1.52 


0.27 


SEM 


(12) 


30 


0.26 


0.05 


1.54 


0.28 


ISA 


(13) 


5 


0.37 


0.09 


1.62 


0.29 


EET 


(13) 


24 


0.32 


0.06 


1.46 


0.46 


ESE 


8 


0.34 


0.16 


1.33 


0.49 


ITM 


1 


0.39 


0.10 


0.99 


0.12 



Tabla 1 Valores medios y desviaciones de k y yencon trades para 

cada una de las estaciones de la RAM 

En la tabla 1 se presentan las estaciones acelerograficas y el numero de la microzona sfsmica 
correspondiente, el numero de registros integrables de cada estacion, N, y los valores medios y las 
desviaciones de k y y, encontrados. 

El parametro k oscila entre 0.22 y 0.39, con desviaciones para cada estacion del orden de 0.06, lo 
que habla muy bien de la estabilidad de este parametro. El parametro y oscila entre 0.99 y 1.80, con 
desviaciones del orden de 0.50, lo que lo hace un parametro menos estable que el anterior. 

En algunos registros, aunque se considero que el proceso de integracion cumplio con los 
requisitos establecidos en el procedimiento para procesamiento de los registros, el periodo T 2 , donde 
se presenta el desplazamiento espectral maximo, esta considerablemente alejado del que se reporta 
en la literatura como periodo fundamental de vibracion del deposito de suelo en el sitio. Lo anterior 
significa que estos registros presentan contribucion considerable de aceleracion en periodos altos. 
La presencia de aceleraciones en estos periodos puede provenir de efectos regionales, que se 
descartan porque solo unos cuantos registros lo muestran, y deberian mostrarlo, si no todos los 
registros, si la mayoria de ellos. Una explication mas crefble es que estas componentes de 
aceleracion en periodos altos provienen directamente de las fuentes sismogenicas. 

Una mirada mas detallada a estos registros, que son los mismos en todas las estaciones, 
coinciden en el hecho de que todos provienen de fuentes profundas, especialmente de la zona de 
Benioff del Viejo Caldas. Esta particularidad hace pensar que efectivamente el fenomeno descrito 
es real y debe tenerse en cuenta. 
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Si se considera que el desplazamiento espectral en estos periodos altos alcanza el mismo nivel 
que el desplazamiento espectral en el periodo fundamental del sitio, hipotesis razonablemente 
conservadora, y ademas, se consideran los porcentajes de participacion de la zona de Benioff del 
Viejo Caldas en la amenaza sismica de Medellin para las intensidades de diseno y control de danos. 
es posible estimar el periodo de las ondas hasta las cuales se cubre con cada parametro y. 

En los resultados que se presentan a continuacion, se considera que el periodo fundamental del 
sitio coincide con el periodo T 2 que especifica la microzonificacion sismica de Medellin para cada 
una de las zonas. Igualmente, se considera, basados en el estudio de amenaza sismica del proyecto 
de microzonificacion sismica de Medellin, que los porcentajes de participacion de la zona de 
Benioff del Viejo Caldas en la amenaza sismica de Medellin son de 33% y 40% para las 
intensidades de control de danos y de diseno, respectivamente. Finalmente, se considera en este 
calculo un parametro k igual a 0.30, que aunque no es el menor reportado en la tabla 1. si puede 
considerarse como su limite inferior, que es el que hace mas desfavorable este calculo. Mas 
adelante se concreta esta propuesta. 





T cubierto 


ZONA 


Sa m „ 


t 2 


Sd m , x 


y=1.50 


Y=1.00 


y=0.50 


1 


0.70 


0.60 


6.26 


2.20 


4.20 


29.40 


2 


0.80 


0.40 


3.18 


1.46 


2.80 


19.60 


3 


0.80 


0.70 


9.74 


2.56 


4.90 


34.30 


4 


0.50 


0.65 


5.25 


2.38 


4.55 


31.85 


5 


0.60 


0.60 


5.37 


2.20 


4.20 


29.40 


6 


0.50 


0.50 


3.11 


1.83 


3.50 


24.50 


7 


0.60 


0.60 


5.37 


2.20 


4.20 


29.40 


8 


0.55 


0.75 


7.69 


2.74 


5.25 


36.75 


9 


0.70 


0.55 


5.26 


2.01 


3.85 


26.95 


10 


0.80 


0.50 


4.97 


1.83 


3.50 


24.50 


11 


0.75 


0.65 


7.87 


2.38 


4.55 


31.85 


12 


0.80 


0.70 


9.74 


2.56 


4.90 


34.30 


13 


0.80 


0.50 


4.97 


1.83 


3.50 


24.50 


14 


0.60 


0.55 


4.51 


2.01 


3.85 


26.95 



Tabla 2 Periodos cubiertos para la intensidad de diseno con diferentes 

parametros y, si se presentan componentes de periodos altos 
provenientes de la zona de Benioff del Viejo Caldas que 
generen desplazamientos espectrales iguales a los que se 
presentan en el periodo T 2 de la zona. 





ESPECTROS DE DISENO SISMICO PARA LA CIUDAD DE MEDELLIN 



11 





T cubierto 


ZONA 


Sa m „ 


t 2 




y=1.50 


Y=1.00 


Y=0.50 


1 


0.70 


0.60 


6.26 


4.90 


14.00 


326.67 


2 


0.80 


0.40 


3.18 


3.27 


9.33 


217.78 


3 


0.80 


0.70 


9.74 


5.72 


16.33 


381.11 


4 


0.50 


0.65 


5.25 


5.31 


15.17 


353.89 


5 


0.60 


0.60 


5.37 


4.90 


14.00 


326.67 


6 


0.50 


0.50 


3.11 


4.08 


11.67 


272.22 


7 


0.60 


0.60 


5.37 


4.90 


14.00 


326.67 


8 


0.55 


0.75 


7.69 


6.12 


17.50 


408.33 


9 


0.70 


0.55 


5.26 


4.49 


12.83 


299.44 


10 


0.80 


0.50 


4.97 


4.08 


11.67 


272.22 


11 


0.75 


0.65 


7.87 


5.31 


15.17 


353.89 


12 


0.80 


0.70 


9.74 


5.72 


16.33 


381.11 


13 


0.80 


0.50 


4.97 


4.08 


11.67 


272.22 


14 


0.60 


0.55 


4.51 


4.49 


12.83 


299.44 



Tabla 3 Periodos cubiertos para la intensidad de control de danos 

con diferentes parametros y, si se presentan componentes 
de periodos altos provenientes de la zona de Benioff del Viejo 
Caldas que generen desplazamientos espectrales iguales a los 
que se presentan en el periodo T 2 de la zona. 



El calculo anterior es valido solo si el valor del parametro k es menor que el porcentaje de 
participacion de la zona de Benioff del Viejo Caldas en la amenaza sfsmica de Medellin. De las 
tablas 2 y 3 es claro que para valores de ymenores o iguales a 1.0 ya se tiene un cubrimiento de 
estas ondas hasta periodos suficientemente altos. 

Para tratar de cubrir lo que pasari'a en el caso de registros de altas intensidades, con los que no se 
cuenta todavfa hoy, se considera que al menos el parametro k puede estar relacionado con la forma 
de los espectros de aceleracion en la zona de periodos bajos e intermedios, de la que por experiencia 
debido al trabajo con espectros de aceleracion se tiene al menos buena sensibilidad del parametro 
que controla esta zona del espectro, Sa mw JA max . En lo que sigue se aclara esta propuesta y se define 
un criterio mas robusto que permita proponer el parametro k. 

Aunque los resultados arrojados por el procesamiento de los acelerogramas registrados por la 
RAM no muestran una correlacion clara y definitiva entre el cociente Sa mr JA mru y el parametro 1 Ik 
estimado para cada una de las zonas, el conocimiento que se tiene del fenomeno indica que deberfa 
haberlo. El cociente Sa m JA mcu indica que tan concentrado o que tan piano es el contenido 
frecuencial de la serial que se representa: en zonas de suelos blandos este cociente puede llegar a 
tener valores del orden de 7.0, indicando una gran concentracion de las frecuencias de la serial 
alrededor de la frecuencia fundamental del sitio. En suelos duros este factor esta alrededor del 
famoso 2.5, que durante mucho tiempo sirvio de base para definir la forma de los espectros de 
diseno alrededor del mundo. El inverso del parametro k representa un cociente semejante al 
discutido en el parrafo anterior, pero esta vez en el caso del espectro de desplazamiento, 1 Ik - 
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Sd ma JD max , que tambien tiene relacion con la concentration del contenido frecuencial de la senal 
que se quiere representar. 

Una explication para la falta de evidencia, en los datos disponibles en Medellin, de la 
correlacion entre 1 Ik y el cociente Sa ma JA max que se acaba de discutir, es la escasez de registros de 
contenido frecuencial relativamente piano. Por razones de todos conocidas, los sismos en estas 
zonas duras sistematicamente registran aceleraciones maximas del suelo, Amax, mucho menores 
que las que se registran en las zonas de suelos mas blandos, lo que hace que el numero de registros 
integrables en estas zonas duras sea muy escaso en comparacion con los que se tienen en las zonas 
de suelos blandos. 

Aun sin evidencia experimental, se considera pertinente tener en cuenta la correlacion 
mencionada. Se propone, entonces, que haya una relacion directa entre el cociente Sa m JA max y 1 Ik, 
de tal manera que para valores del primero menores o iguales a 2.5 se considere un valor para 1 Ik 
igual a uno, es decir, un espectro de desplazamiento piano a partir del periodo T 2 . Para valores del 
cociente Sa ma JA miu mayores a 2.5 y hasta 4.6, que es el maximo valor reportado por la 
microzonificacion sfsmica de Medellin en los espectros de control de danos, se considere un valor 
de 1 Ik creciente a partir de uno y hasta el valor lfmite razonable reportado por los datos, 1/0.30. 

La consideration anterior elimina en muchas zonas la posibilidad de decidir sobre el valor de y 
para cubnr contenidos de aceleracion en periodos altos provenientes de la zona de Benioff del Viejo 
Caldas. En particular, se elimina esta posibilidad en todas aquellas zonas con valores de k mayores 
a 0.4, para el caso de los espectros de diseno, y en zonas con k mayores a 0.33, en el caso de los 
espectros para control de danos. 

Lina ultima verificacion a los espectros por proponer se hace atendiendo la observation de 
algunos autores que sugieren que la rama descendente de los espectros de aceleracion debe 
comenzar con una cafda proporcional a 1 IT, para representar la zona de velocidad constante, 
terminada la cual hay una zona de transition a una cafda proporcional a 1/7 2 , que fmalmente 
representa la zona de desplazamiento constante. Los mismos autores que sugieren estas formas 
espectrales justifican los espectros actuales con cafdas proporcionales a 1/7’ al hecho de que la zona 
de velocidad constante se ha ido ‘ alargando” hasta el punto que algunos reglamentos, como la 
norma NSR-98, la proponen como la rama final de los espectros de aceleracion. En estos espectros, 
es claro que el desplazamiento espectral maximo se presenta al final de la zona de velocidad 
constante y no al final de la zona de aceleracion constante como se sugiere en esta propuesta. 

La forma espectral descrita en el parrafo anterior es caracterfstica de suelos duros, que en 
presencia de estratos superflciales de suelos mas blandos se ve completamente modificada por la 
luncion de transferencia, y es esta mezcla la que fmalmente domina el panorama en muchas zonas 
uibanas del mundo y en particular en la ciudad de Medellin. No obstante lo anterior, se hara la 
verificacion que se describe a continuation. 

En lo que sigue se considera que la zona de velocidad constante del espectro en el basamento 
rocoso se extiende hasta un periodo de 1.2 seg,. como lo sugieren los espectros de aceleracion de la 
estacion en roca para Medellin. Para conseguir la aceleracion en este periodo estructural pero a 
mvel de la superficie del suelo, la aceleracion en este periodo a nivel de la roca se debe multiplicar 
por la funcion de transferencia sobre espectros de respuesta de aceleracion, que es funcion del 
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cociente entre el periodo sobre el que se quiere la funcion de transference (el de las ondas 
incidentes), en este caso 1.2s, y el periodo fundamental del deposito de • suelo, y del 
amortiguamiento del suelo. 

Debido a que en el caso de Medellin el deposito de suelo con periodo mas largo es del orden de 
0.60 seg, como lo sugieren las tablas 4 y 5, lo que resulta en un cociente entre el periodo 
fundamental del suelo y el periodo de las ondas incidentes, 0.60/1.2, igual a 0.50, bastante alejado 
de 1.0, donde se presentarian las amplificaciones maximas, y por las razones expuestas en parrafos 
anteriores y que tienen que ver con la forma de la funcion de transferencia, se considera 
conservadoramente un factor de amplificacion sobre el espectro de respuesta de aceleraciones de 
1.2 para el cociente Sa max /A max igual a 4.6, e igual a 1.0 para el cociente Sa mw JA max igual o menor a 
2.5. Para valores intermedios del cociente Sa ma JA max se hace interpolation lineal. Ahora, debido a 
la cercam'a del cociente entre periodos al valor 0.0, donde la funcion de transferencia tiende a 1 .0 
independientemente del amortiguamiento, el factor de amplificacion se calcula con la misma 
formula tanto para los espectros de control de danos como de diseno. 

La aceleracion espectral a nivel del basamento rocoso y correspondiente a un periodo estructural 
de 1.2s se calcula suponiendo que la aceleracion espectral maxima, Sa mnx , es igual a la aceleracion 
maxima a nivel de roca, A max , multiplicada por el factor famoso de 2.5. Ademas, se supone que a 
nivel de roca, el espectro cae proporcional a 1 IT (zona de velocidad constante) entre el periodo T 2 y 
1.2 s. Para el espectro de control de danos a nivel de roca se supone que el periodo de quiebre T 2 es 
igual a 0.3s, y para el espectro de diseno igual a 0.4s. En estas condiciones, la aceleracion espectral 
en 1.2s a nivel del basamento rocoso resulta igual a 0.019g para el espectro de control de danos, e 
igual a 0.1 25g para el espectro de diseno. 

Basados en el mini mo valor que para y resulta de los datos disponibles en cada una de las 
microzonas sismicas (ver tabla 1), en la restriction que en unas pocas zonas introduce la necesidad 
de considerar las ondas de periodos altos provenientes de la zona de Benioff del Viejo Caldas, y 
finalmente, en la condition de que la aceleracion espectral en 1.2s sea igual o superior a la 
aceleracion espectral para este periodo calculada con el procedimiento descrito en el parrafo 
anterior, se llega a los parametros que se muestran en las tablas 4 y 5 que defmen los espectros de 
control de danos y diseno, respectivamente, en las diferentes zonas definidas por la 
microzonificacion slsmica de Medellin. 

Es importante observar que para poder cumplir el ultimo requisito propuesto es necesario alargar 
ligeramente, 0.1s, los periodos de quiebre T 2 de los espectros de control de danos y de diseno en las 
microzonas 2 y 6. En la zona 2 el problema se presenta porque los periodos de quiebre T 2 originales 
son de solo 0.3s para control de danos y de 0.4s para diseno, lo que hace que los espectros 
propuestos caigan fuertemente entre estos periodos y 1.2s. En la zona 6 el problema no es solo que 
los periodos T 2 sean cortos sino que la meseta de aceleracion constante esta muy baja, lo que 
produce aceleraciones espectrales muy bajas para 1.2s. 

A manera de comparacion, se incluyen a continuation de las tablas, las graficas de los espectros 
definidos en la microzonificacion sismica de Medellin, en linea punteada, y los definidos por este 
trabajo, en linea continua, para cada una de las microzonas. 
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Era 




Samax 


Samax/Amax 


T1 


T2 


k 


Y 


i 


0.05 


0.23 


4.60 


0.10 


0.50 


0.30 


1.00 


2 


0.08 


0.22 


2.75 


0.10 


0.40 


0.78 


1.00 


3 


0.07 


0.25 


3.57 


0.10 


0.50 


0.46 


1.00 


4 


0.05 


0.18 


3.60 


0.10 


0.60 


0.45 


1.00 


5 


0.06 


0.22 


3.67 


0.10 


0.50 


0.44 


1.00 


6 


0.05 


0.14 


2.80 


0.10 


0.50 


0.75 


1.00 


7 


0.06 


0.22 


3.67 


0.10 


0.50 


0.44 


1.00 


8 


0.08 


0.18 


2.25 


0.10 


0.65 


1.00 


1.00 


9 


0.06 


0.23 


3.83 


0.10 


0.40 


0.40 


0.30 


10 


0.09 


0.25 


2.78 


0.10 


0.40 


0.76 


1.00 


11 


0.06 


0.23 


3.83 


0.10 


0.50 


0.40 


1.00 


12 


0.06 


0.25 


4.17 


0.10 


0.65 


0.35 


1.00 


13 


0.06 


0.25 


4.17 


0.10 


0.40 


0.35 


0.30 


14 


0.05 


0.14 


2.80 


0.10 


0.50 


0.75 


1.00 



Tabla 4 Parametros que definen los espectros de control de danos 

para las microzonas sfsmicas de Medellin 



ZONA 


Amax 


Samax 


Samax/Amax 


T1 


T2 


K 


y 


1 


MEM 






0.10 


0.60 


0.91 


1.00 


2 


H 




2.35 


!■ 


■ 


1.00 


1.00 


3 


0.30 


0.80 


2.67 


0.20 


0.70 


0.84 


1.00 


4 


0.23 


0.50 


2.17 


0.10 


0.65 


1.00 


1.00 


5 


0.18 


0.60 


3.33 


0.10 


0.60 


0.52 


0.30 


6 


0.18 


0.50 


2.78 


0.10 


0.60 


0.76 


0.00 


7 


0.18 


0.60 


3.33 


0.10 


0.60 


0.52 


0.30 


8 


0.23 


0.55 


2.39 


0.10 


0.75 


1.00 


1.00 


9 


0.26 


0.70 


2.69 


0.10 


0.55 


0.83 


1.00 


10 


0.38 


0.80 


2.11 


0.10 


0.50 


1.00 


1.00 


11 


0.26 


0.75 


2.88 


0.10 


0.65 


0.70 


1.00 


12 


0.26 


0.80 


3.08 


0.15 


0.70 


0.61 


1.00 


13 


0.26 


0.80 


3.08 


0.10 


0.50 


0.61 


0.10 


14 


0.20 


0.60 


3.00 


0.10 


0.55 


0.64 


0.00 



Tabla 5 



Parametros que definen los espectros de diseno 
para las microzonas sfsmicas de Medellin 
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Figura 3 



Espectros de diseno y de control de dano propuestos 

en la microzonificacion si'smica de Medellin y en este trabajo. 
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Figura 3 Espectros de diseno y de control de dano propuestos 

en la microzonificacion si'smica de Medellin y en este trabajo. 

Los espectros de aceleracion para diseno y control de danos propuestos en este trabajo reducen 
sustancialmente las aceleraciones a partir del periodo de quiebre 7 2 , y como es de esperarse, muy 
rapidamente son inferiores a la aceleracion espectral minima propuesta para diseno en la 
microzonificacion si’smica de Medellin, A nu J2. Por considerar que esta limitacion en las 
aceleraciones exigidas para diseno es razonable, se propone entonces que aunque la aceleracion 
minima para diseno mantenga el li'mite establecido, el calculo de los desplazamientos para efectos 
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de control de derivas se haga con base en el espectro propuesto de aceleraciones, o en el de 
desplazamientos derivado de este ultimo. 



CONCLUSIONES 



Es claro que los espectros de aceleracion que actualmente propone la microzonificacion slsmica 
de Medellin violan en la zona de periodos altos los principios de la Dinamica Estructural: los 
espectros de desplazamiento resultantes crecen indefinidamente, lo que resulta en un 
desplazamiento maximo del terreno infinito. Este hecho aunado a los llmites que actualmente 
controlan el diseno slsmico de estructuras, como son los llmites de derivas, resultan en sobre 
disenos que superan en la mayorla de los casos cualquier llmite racional que se pueda fijar a un 
factor de seguridad. 

En este trabajo, basados en los principios de la Dinamica Estructural y en los registros hasta 
ahora obtenidos por la Red Acelerografica de Medellin, desafortunadamente de muy bajas 
intensidades todos ellos (2 cm/seg 2 <A mcLX <\() cm/seg 2 ), se proponen espectros para diseno y control 
de danos en cada una de las zonas defmidas por la Microzonificacion Slsmica de Medellin. En esta 
propuesta no se involucran las partes de los espectros en las zonas de periodos bajos o ascendente 
del espectro de aceleraciones, ni la zona de periodos intermedios o de amplificaciones maximas. 
Mas aun, se consideran estos valores como ciertos, y se usan para deftnir a partir de estos la forma 
de calda de los espectros. 

Como se menciona en la introduccion, a medida que se recabe mas information acelerografica 
se puede acudir a metodologlas mas acordes con lo que se busca definir, como puede ser el ajuste de 
leyes de atenuacion, preferiblemente fundamentadas en principios de transmision de ondas, y a 
partir de estas ajustar las formas de espectros de amenaza uniforme en la zona de periodos altos. 
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Resumen 

En la actualidad existen muchos puentes, construidos hace veinte o treinta anos, que muestran reducciones 
diversas en su capacidad resistente. Para estas estructuras son necesarias herramientas de evaluacion 
simplificadas y confiables de su riesgo, que permitan obtener una estimacion de su futuro comportamiento. En 
busca de ese objetivo se propone un metodo de analisis a partir de un modelo dinamico y de calculo de la 
maxima inercia degradada de las pilas del puente. Utilizando este metodo se evalua la vulnerabilidad sismica de 
un puente de hormigon armado existente, construido en los anos setenta bajo requisitos de diseno sismico ahora 
obsoletos. Para esta estructura se definio un posible escenario sismico, para el cual se generaron familias de 
registros artificiales. Ademas, se realizo una simulacion por el metodo de Monte Carlo, considerando las 
incertidumbres asociadas a la peligrosidad sismica y a los parametros estructurales. A traves de los resultados 
obtenidos se comprueba la bondad del proceso simplificado propuesto. 



Abstract 

Actually, many of the existing bridges, built twenty or thirty years ago, show different reductions in their 
resistant capacity. For these structures, simplified and reliable evaluation tools that could be applied in the 
evaluation of their future behavior are necessary. Therefore, a method of analysis based on a simplified dynamic 
model and in the calculation of the maximum damaged inertia of the bridge piers, is proposed. Using this 
method the evaluation of the seismic vulnerability of an RC existing bridge, built under now obsolete seismic 
design requirements, is proposed in the paper. For the structure site a possible seismic scenario was defined, for 
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which families of artificial record were generated. Furthermore, to evaluate the vulnerability of the structure a 
Monte Carlo simulation was performed, considering the uncertainties in the seismic hazard and in the structural 
parameters. The obtained results allow assessing the goodness and reliability of the proposed simplified process. 



1. INTRODUCTION 

Un gran numero de los puentes de los actuales sistemas de transporte fue construido hace veinte o 
treinta anos. Algunas de estas estructuras presentan ahora diferentes grados de deterioro, que es 
necesario evaluar de manera adecuada para tomar decisiones para su futuro refuerzo y 
rehabilitacion (1) . 

En la estimacion del comportamiento actual de estructuras existentes es ya comun la 
determinacion de su vulnerabilidad sfsmica mediante curvas de fragilidad o matrices de probabilidad 
de dano. En estudios probabilistas como estos, los analisis multiples, que consideran las 
incertidumbres inherentes al problema, requieren el uso de metodos de calculo confiables, simples y 
economicos. Asf, desde hace mas de una decada, se han propuesto numerosos metodos de estimacion 
del dano, en su mayorfa para la caracterizacion del comportamiento de edificios <2) . Actualmente, 
existen algunos metodos de evaluacion para puentes, englobados en cuatro grupos principales a 4) : (1) 
determinacion del fndice de vulnerabilidad, (2) utilization de modelos estructurales de evaluacion, (3) 
estimaciones basadas en el juicio de expertos y (4) defmiciones a partir del analisis estadfstico de 
datos reales. El primer metodo esta basado en estimaciones simples (muchas veces visuales) que 
buscan una primera clasificacion de las estructuras con mayores problemas. Por su parte, los metodos 
basados en modelos estructurales suelen ser herramientas que proporcionan gran cantidad de 
resultados, pero cuya confiabilidad depende de su calidad y facilidad de aplicacion. Finalmente, las 
estimaciones basadas en el juicio de expertos requieren disponer de un buen numero de profesionales 
con conocimientos y experiencia comprobada, mientras que las estimaciones estadisticas de datos 
reales solo se aplican en zonas de sismicidad moderada o severa con information suficiente. 

Con el objetivo de proporcionar una herramienta de analisis simplificada y confiable para 
la estimacion de la vulnerabilidad de puentes de hormigon armado de eje simple de pilas, mas 
alia de la simple determinacion del fndice de vulnerabilidad y sin llegar a analisis demasiado costosos 
con una gran cantidad de resultados, se propone una herramienta de caracterizacion del dano (4) . En 
este artfculo se describe primeramente el procedimiento de analisis propuesto. En este, la respuesta 
dinamica no lineal de la estructura se obtiene a traves de una formulacion incremental (paso a paso). 
un procedimiento de solution iterativa hasta una convergencia adecuada, de acuerdo con la tecnica de 
Newton-Raphson, y un algoritmo de integration en el tiempo englobado en el metodo de Newmark. 
En el proceso de linealizacion de la ecuacion de movimiento aplicado en este tipo de analisis se 
actualiza en cada paso la matriz de rigidez tangente de la estructura. Con esto, la formulacion lineal 
simplificada propuesta se utiliza no solo en la calibration del modelo, sino tambien en la 
determinacion del predictor de carga externa en cada iteration e incremento de tiempo. 

La metodologfa simplificada propuesta es utilizada como herramienta de analisis en la estimacion 
de la vulnerabilidad sfsmica del puente de hormigon armado Warth (perteneciente a la tipologfa de 
estudio), situado sobre la autopista Viena-Graz, en Austria. Conforme a las normas de diseno vigentes 
en la zona de ubicacion, el puente requiere estudios que indiquen su estado actual para la toma de 
futuras decisiones. 

Para conocer la vulnerabilidad del puente Warth se realizo una simulacion por el metodo de Monte 
Carlo, en la que se consideraron las incertidumbres inherentes a la excitation sfsmica y a las 
propiedades estructurales. A traves de esta simulacion se determinaron curvas de fragilidad de la 
estructura y de sus pilas, con lo que se define el posible dano de la estructura. 

Aunque existe una variada gama de metodos de evaluacion de la respuesta no lineal de puentes, la 
metodologfa propuesta representa un camino simple, confiable y de bajo costo computacional. Por lo 





